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Abstract 

We study the fermion propagator in a spatially varying classical background field, and show 
that, contrary to common wisdom, it may get nontrivial gradient corrections already at the 
first order in derivative expansion. This occurs whenever the fermion self-energy acquires a 
spatially (or temporally) varying pseudoscalar term, a simple example of which is given by 
a complex mass term rh{x) = mji + i'~fr,mj. Such effective mass terms arise for example in 
extensions of the Standard Model during the electroweak transition, and they are crucial in 
providing the CP-violation necessary for electroweak baryogenesis. 



1 Introduction 

Often in applications of quantum field theories, it is desirable to treat a part of the problem 
semiclassically and compute the quantum corrections in a controlled expansion around the 
classical solutions. This applies to a broad spectrum of problems, including for example the 
Fermi liquid theory, as well as the study of the expanding and cooling plasma in the heavy 
ion collisions 0] . In these, and in many other cases, it is often sufficient to assume that 
the background is very slowly varying and to treat the problem at leading order in gradients. 
Such an adiabatic (quasiparticle) approximation may make sense if the physical effect of 
the background one is interested in is already present at the trivial order. However, there 
are situations where a semiclassical approximation is justified, but the system displays an 
important physical effect only at a nontrivial order in gradients. An important example of 
this occurs at the electroweak transition, where the problem is to study particle propagation 
in the presence of a spatially varying, CP-violating Higgs condensate. The CP-violating 
effects of the condensate on the plasma dynamics arise only at the first nontrivial order in 
gradient expansion p, ^, ^, ^. The induced perturbations in the particle and antiparticle 
densities, coupled with the weak anomaly, lead to a local baryon production mechanism, 
capable of creating the observed matter-antimatter asymmetry of the Universe [0, ^. 

Unfortunately, no consistent derivation of the semiclassical quantum Boltzman equation 
for a fermionic field beyond the trivial order in gradients exist in the literature. Moreover, 
there exists a well known no-go theorem according to which first nontrivial corrections 
to fermion propagator arise only at second order in derivative expansion. This is in apparent 
contrast with the classical force mechanism (CFM) for baryogenesis according to which 

the nontrivial effects on fermion propagation are felt already at the first order in gradients. In 
this letter we show that the no-go theorem does not hold universally, but instead, whenever 
the fermion self energy contains a pseudoscalar term, there are corrections to the propagator 
already at the first order in gradients. A simple example of such a term is provided by a 
Dirac field with a complex spatially (or temporally) varying mass term of the form 



This term could be of course replaced by generic scalar and pseudoscalar operators; we retain 
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the present notation however, because of its clear physical implications. Namely, this type 
of term has precisely the form relevant for baryogenesis, where such an effective mass term is 
induced by a coupling of the fermion with a spatially varying complex scalar field condensate 
during a first order phase transition. Here we shall consider a simple toy model with the 
lagrangian 

£ = i^^^ - ^^mipR- i/jR m* ipL + 'Cint, (2) 

where m = tur + imj is a complex mass. The detailed form of the interaction lagrangian /^int 
is not relevant for us here. Realistic cases studied in literature in connection with the CFM 
baryogenesis include the two doublet extension of the standard model and the chargino 
sector of the minimal supersymmetric standard model 0. 

This letter is a part of a series P, ^ |10[ dedicated to derivation of semiclassical transport 
equations for quantum fields in varying backgrounds beyond the trivial order in gradients. 
We first explicitly construct the fermion propagator for the theory (0) to the first nontrivial 
order in gradients. Corrections appear as higher order poles multiplied by derivatives of the 
complex phase of the mass. The position of the poles is not shifted however, and the phys- 
ical content of the correction terms only becomes apparent after consideration of spectral 
integrals over test functions by a technique developed in . It will be shown that the spec- 
tral function, while not expressible as a sum of ordinary delta-functions even in the on-shell 
limit, projects a test function to a sharp, but shifted energy shell. We then show that this 
shell corresponds to the physical semiclassical dispersion relation of a standard WKB wave 
function, and describe how the current computed using the spectral function representation 
can be related to motion of WKB wave packets. 



2 Propagator 

Our first task is to solve the retarded and advanced propagators for the fermionic field ip 
in the weak coupling approximation and in an expansion in gradients of the slowly varying 
mass terms. We will work in the Keldysh closed time contour (CTC) formalism 



2 



and the basic quantity of interest is the path-ordered propagator 

iGt{x,y) = {Tcmy)^l,{x)]). 



(3) 



Writing the Schwinger- Dyson equations for Gc in the real time formulation implies the 
following formally exact equations for the retarded and advanced propagators G^''^ in the 
Wigner representation 

e-'''{{iG,)-'}{iG^^\k-X)} = l, (4) 

where ^{f}{g} = \{,dxf ■ dkQ — dkfdxg) and the tree level propagator is given by 

itG'-^T' = rriRiX) - t^.mjiX) - T.'''^ ± is^e, (5) 

where s^^ = sign(ci;), k denotes the canonical 4-momentum and X = {x + y)/2 the average 
position. From now on we shall assume a weak coupling limit and neglect the self-energy 
Y;r,a _ xi^zpituF; in particular we thus consider the on-shell limit F— >0. We will also suppress 
the ie-prescription which differentiates between C'" even when F^O, which is introduced 
as a mnemonic tool to maintain information about the integration path required to fully 
specify how G acts operationally. A formal solution to equation (P) can easily be found in 
the gradient approximation. Indeed, to first order in gradients we have: 

iG = tGo - tGo tO{itGo)-'}tGo{itGo)-'} tGo, (6) 

where the free propagator is given by 

'Go = ^, + mR{X) - zj^miiX)). (7) 

/e^ — \m\ 

After some algebra one can write the expression into the following covariant form: 

iG = iGo + ] , [tik ■ dx - J^^^)\m\e-''^' + tVI'^x^] • (§) 

There is clearly a nonvanishing correction to the propagator at the leading nontrivial order. 
This is not in disagreement with the no-go theorem by because, as we shall see, all 
nontrivial consequences of (|]) are proportional to the complex pseudoscalar mass (or more 



generally to a pseudoscalar self-energy function), which was not considered in Ref. An 
important point to observe about (H) is that the pole of the propagator is not shifted, but 
instead a new second order pole, proportional to gradients of m has appeared. This behaviour 
is analogous to one previously observed in the case of a scalar field 

For simplicity we shall from now on restrict ourselves to the case of a planar symmetry 
with only spatially varying mass: m = mlx^). In this case one can Lorentz boost to the 

■II 

simplifies to 
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frame in which kn =0, so that u ^ u = s^^^u^ — kj?)^^^. With these definitions Eq. 
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ic:;(m'^7%5 - im^7°) + \mfe' S'\ (9) 
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where we combined various 7-matrices into the spin operator = 7*^7^7^, which measures 
the spin in the ks-direction and in the frame in which ky = 0, and we used the definition z = 
k'^ — \mf. Allowing z a small complex value, one sees that yields G"" {G°') for lm{z) < (> 
0). A further very important simplification arises from the fact that pseudoscalar interactions 
conserve spin (one can show that the spin-operator commutes with the propagator G). 
The problem then becomes diagonal in spin and, by replacement S^^s, the 4x4 structure 
of d^) becomes block-diagonal, so that the problem can be reduced to a 2 x 2 problem. In 
the 2x2 formalism, it is more convenient to use the Pauli matrices for the spinor algebra, 
which can be effected by the following identifications: 

where we have taken the chiral representation for the 7 matrices. In this way (|^) can be 
recast in a particularly simple form: 

z 

— — ulm'ja^ — m'^a'^) + Iml^^^' 
z^ L J 

= 6o + b-(T. (11) 



From the last expression in particular it is evident that there exists a basis in which the prop- 
agator G is diagonal Gs^Gd,s = + s sign(tD)|b|cr^. In the Dirac notation the propagator 
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then reads 

i'yGd = s^\— ^ \ [1 - Sc;,S j , (12) 

where = sign(cc;). Physically the equation ([T^) means that the states with equal ss^j prop- 
agate identically. Formally, this is expressed by the projector (1 — S(2,S'^7^)/2, which identifies 
the states with equal entries on the diagonal. In particular this means that particles and 
antiparticles with opposite spin have identical spectral functions. However, for particles and 
antiparticles of same spin (helicity) the spectral functions differ by an amount proportional 
to the CP-violating angle 6'. This is how in the context of baryogenesis CP-violating back- 
grounds induce bias in the particle and antiparticle distributions. In order to understand 
how this bias propagates in a plasma, it is then necessary to study the relevant quantum 



trasport equations with the effect of collisions included |T3 

3 Spectral integrals 

To study the physical consequences of the propagator (||), we now consider spectral integrals 
over test functions |p , representing some generic observables, such as the generalized particle 
distribution function. We define: 

G[T] = Tr — / duj^'^AT, (13) 



271 Jo 

where A = |(G"" — G") is the spectral function, T is a test function, and the trace is 
taken over the Dirac indices. Because scalar and pseudoscalar backgrounds considered here 
conserve spin, we can take T to be diagonal in the spin space and a scalar function in a spin 
2x2 block. For backgrounds that violate spin this analysis would have to be generalized 
to include spin mixing. The 7°-factor in Eq. (0) is added to make the projection operator 
explicitly hermitean. The spectral function A carries information on the physical spectrum 
of excitations in the system. In particular to zeroth order in gradients it becomes a simple 
on-shell projector (c/. Eq. (|^)): 

A^ [d{uj-iUo)-d{iu + uJo)\,. (14) 



^Note that we are assuming the on-sheh hmit and neglecting possible thermal and/or vacuum corrections, 
so that Y."^''' = 0. 
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where uq = (A;| + |mp)^/^. In a spatially varying background the spectral function can no 
longer be expressed in terms of 5-functions projecting onto the physical energy shells. As 
discussed in [§, it can however be written as a sum of projectors onto complex shells, but 
this is not the approach we will pursue here. Instead, by extending the momentum to the 
complex plane, the integral (13) can be turned into a contour integral over the variable 

G[T] = Tr / dz^^l^, (15) 

47r JCn . /,.>2 I ^ 



'Co ^u^ + z 

where G is the propagator (H), and Co is a contour encircling the pole ai z = {uj = uq 
once clockwise. We can then compute (|T^) by making use of the residue theorem: 



G[T] = Tr Res,=o ^' ^ ' . (16) 

Because the highest pole in the expression (§) is of the second order, only the terms up to 
the linear order in z in the expansions of various terms appearing in ([TBD contribute. In 
particular the gradient expansion in dk^T then terminates at first order. That is to say that 
all 9[^T-terms for / > 1 are of at least second order in gradients. An important consequence 
of this truncation is that a complete description of a fermionic plasma, consistent to first 
order in gradients, requires at most two independent distribution functions, in contrast to 
the four distribution functions required for a full description of the bosonic case at the lowest 
nontrivial order in gradients |l^ . 



It is now a simple matter to compute the residue in Eq. (|T6|) of the propagator (|TT]) or 
(p!^). (For notational simplicity from now on we shall denote u hj u.) To the order we are 
working 

GIT] = %+'-^^^^{%-uoV) 

= (17) 
where Tq = T(co'o), %' = {d^T){u!o), and we have defined the semiclassical energy as 

Use =UJo — ^ — . (18) 

ZuJq 

The test function thus still gets projected to a sharply defined energy shell, which is however 
shifted with respect to the lowest order shell uq. 
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4 Current 



As a physically motivated application of the above, we now consider a current of spin states. 
In the present formalism this can be written as follows: 

jTW = Tr/ ^^^YG<{k■X)P,, (19) 

where = (1 + sS'^)/2 is the spin projector and G^{k; X) is the quantum Wigner function. 
Assuming the decomposition iG'^{k] X)Ps = APsUs, where are generalized particle distri- 



bution functions of spin states, and inserting this expression into (|l^) it is a simple matter 
to show that 

d^k 



(27r)3^, 

jt(X) = 2/-f^^/°(k;X), (20) 

where the (momentum space) semiclassical distribution function /|"^(k;X) = ^^(co'sc, k; X) 
is the projection of Ug onto the semiclassical shell u = Usc, while the distribution function 
(k; X) = ns{uJo, k; X) is the projection of onto the unshifted shell uj = ujq. Further, in 
Eq. (|20|) we have reintroduced the dependencies on the trivial perpendicular directions ky. 
Finally, changing the variables k^^k^c, uJsq—^uj {dk^ — * duj{dujsc/ dk^)~^) in the first integral, 
and k^-^ko = (cu^ — |mp)^/^, luq-^uj {dk^, dLULUo/k^) in the second integral, the current 



reduces to 



dujd^kw I ijj 



where the distribution functions ff' = fs'^{uj, ky; X) = ns{uj, ksc, ky; X) and /f = /f (cu, ky; X) 
= ns(co', fco, k||; X) are now defined on the energy phase space (ci;,k||;X). The semiclassical 
momentum ksc = ^sc(^; k||; X) in Eq. (^) is the one obtained by inverting the semiclassical 
energy-momentum relation ([T8|) : 

^-^^^^)- (-) 

Of course, in equilibrium the spatial currents cancel. However, if the preferred coordinate 
frame is boosted with respect to the plasma rest frame, as is the case for a propagating phase 
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transition front, then there are net currents flowing through the space. One can then deter- 
mine how these induced currents are affected by the presence of a background mass term (|l|). 
Since particles and antiparticles of a given spin (helicity) s and energy Ic^l have different effec- 
tive momenta ksc, the corresponding particle and antiparticle densities are different in the 
region of a spatially varying background. This is how local particle-antiparticle asymmetric 
density and current profiles get created by a spatially varying and CP-violating background, 
which is the phenomenon at the heart of most electroweak baryogenesis mechanisms. 

5 Field theoretical vs WKB approach 

Let us now make a connection between the above results and the dispersion relation obtained 
by using the WKB method. Through a standard calculation ^ 0, one finds 

ssj' , 



fcwkb = ko + -izj—i^ ± -s^o) + /eo = V - \m\ , (23) 

where the arbitrary function a' follows from the reparametrization invariance of the theory 
given by the lagrangian (^ invariant under (global) f/(l)-transformations — > e*"?/^. It is 
then obvious that /c^kb cannot represent a physical quantity. The group velocity 

Vg = (<9^fcwkb)"\ (24) 

corresponding to the stationary phase of the WKB wave packet, is a well defined physical 
quantity however. The physical momentum of a WKB state is then given by /Cphys = ujVg{uj). 
Computing the derivative ( p^ ) one immediately finds the physical reparametrization invari- 
ant WKB dispersion relation 

CSS I Tfi I \ 



which is identical to Eq. (|2^ . With this the semiclassical momentum (|2^) can be reinter- 
preted as the physical momentum of a WKB state. In this sense we have established the 
equivalence between the physical WKB dispersion relation and the field theoretical energy- 
momentum shell corresponding to the propagator (H), or equivalently (0). 

We now return to discuss the current in Eq. (^Tl) . The difference between the WKB 
approach and the field theoretical one presented here is the fact that there are two mutually 
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independent distribution functions /J"^ and /f . In the same spirit as it was done for the 
scalar field in Ref. [10|, we may introduce a coherent quantum density /^*^ = /J"^ — /f, which 
is a dynamical measure of quantum coherence on phase space between the semiclassical and 
the unshifted shell. At the leading order in gradients f^"^ of course vanishes. In analogy to 
the scalar field case we then expect that in the frequent scattering limit f^^ gets suppressed 
and can be consistently neglected. Assuming this is the case, we may then set f^'^ —>■ fs, 
fg fs, so that finally the current (|2T|) reduces to the following form 

This limit is in perfect agreement with the naive WKB result implying the following 
simple physical interpretation: when quasiparticles arrive into a region of a spatially varying 
background, they either speed up or slow down. If they slow down for example, the local 
particle density increasess proportionally to the factor 1/vg. At the same time the flux, 
given by the spatial part of the current, remains unchanged, as it should. To prove ( ^61) 
rigorously requires derivation and solution of the relevant quantum transport equations to 
determine self consistently the distribution functions /l'^ and /°. This problem will be further 
considered in a forthcoming work ||13|. 

While it was to be expected that the field theoretical calculation, unlike the WKB- 
computation, is automatically invariant under the field redefinition ip —>■ e^^'ip, the fact 
that the two methods lead to identical dispersion relations and currents (in the frequent 
scattering limit) is highly nontrivial. Indeed, it should be stressed that the physical WKB- 
shell has emerged from the field theoretical calculation indirectly, not as a simple shift of 
the propagator poles, but rather by the operation of a modified spectral projector on test 
functions, which is effected via the technique of spectral integrals. 



6 Conclusions 

We have studied the fermion propagator coupled with a spatially varying classical background 
field. We represented this coupling by a spatially varying complex mass term, and proved 
that the propagator receives nontrivial corrections already at the first order in the gradients 



9 



of the mass term m. For this result it was crucial that m contains a nonvanishing complex 
pseudoscalar piece. We then explicitly constructed the fermion propagator and the associated 
spectral function to the first order in gradients. The gradient corrections were shown to 
arise, not as shifts in the poles, but rather as additional higher order poles multiplied by 
the gradients of the complex phase of m. We then considered a spectral integral of a simple 
test function, and showed that it nevertheless gets projected to a sharp energy shell, which 
acquires a spin dependent shift with respect to the lowest order shell. The new shell was then 
shown to coincide with the physical dispersion relation of a WKB-state. We also computed 
the current in the field theoretical way, and showed that it can be interpreted in terms of 
moving WKB wave packets. We note that it is straightforward to generalize our results to 
a generic pseudoscalar self-energy term and to temporally varying fields. 

Here we have concentrated on the propagation of fermions in a nontrivial backgrounds. 
The next logical step is to study the transport properties in a plasma containing fermions. 
Drawing on the insights obtained from the study of the scalar case [|lO] , the results presented 



here suggest that a consistent description of transport at nontrivial order in gradient expan- 
sion will require a set of coupled equations for number density and an additional function 
describing quantum coherence in phase space. We also expect that in the frequent scattering 
limit the quantum coherence effects are suppressed, and an effective quasiparticle picture is 
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